The Josephson relation for the superfluid density in the BCS-BEC crossover 
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The Josephson relation for the superfluid density is derived for a Fermi superfluid in the BCS- 
BEC crossover. This identity extends the original Josephson relation for Bose superfluids. It gives a 
simple exact relation between the superfluid density ps and the broken-symmetry Cooper pair order 
parameter Ao in terms of the infrared limit of the pair fluctuation propagator. The same expression 
holds through the entire BCS-BEC crossover, describing the superfluid density of a weak-coupling 
BCS superfluid as well as the superfluid density of a Bose condensate of dimer molecules. 
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I. INTRODUCTION 

Josephson's relationi^ has played an important role 
in developing an understanding of superfluidity since it 
establishes the connection between the two order param- 
eters that are widely used to discuss Bose superfluids: the 
condensate density ric and the superfluid density ps- It 
provides an exact relation between these quantities in a 
Bose superfluid in terms of the infrared behaviour of the 
single-particle Green's function Dn for bosons of mass 
mB- 



= — lim „ „ — r. 

q-o q2i?ii(q,0) 



(1) 



Here, Z?ii(q, 0) is the single-particle Green's function for 
bosons with momentum q at zero Matsubara frequency, 
Vm = 0. 

The simple structure of Eq. ([1]) has enabled a detailed 
analysis to be carried out of the superfluid transition 
in both three— 1^ and two-dimensional'^ Bose superfluids. 
This includes finite-size systems where it has been used 
to study the Berezinskii-Kosterlitz-Thouless transition^ 
in two-dimensional superfluids. In particular, by start- 
ing from a finite two-dimensional system, Ref. |j explic- 
itly showed how a nonzero superfluid density can persist 
in the thermodynamic limit even though the condensate 
density vanishes. 

Recently, Holzmann and Baym have extended the orig- 
inal phenomenological arguments of Josephson and given 
a microscopic proof of Eq. ([1]) using diagrammatic per- 
turbation theoryj^ This proof extends earlier discussions 
by Bogoliubov^ Gavoret and Nozieres;^ and Hohenberg 
and Martin^ at T = where ps = mn. A discussion 
of the Josephson relation for a Bose superfluid at finite 
temperatures is given by Griffini^ using the dielectric di- 
agrammatic formalism (see also Wong and Gould^^). 

Following Holzmann and Baym, in this paper the anal- 
ogous exact relation for a two-component Fermi super- 
fluid is proven. This gives the relationship between the 
superfluid density and the order parameter Aq that rep- 
resents the Bose-condensate of Cooper pairs of fermions. 
The Josephson relation derived in this paper is analyzed 



in the BCS-BEC crossover-i^iHiil picture of Fermi super- 
fluids, widely studied in recent years in the context of 
ultracold atomic gasesi ^^d^d'^d^ In the BEC limit of this 
crossover, where the attractive interaction between the 
two species of fermions is strong, the Cooper pairs re- 
duce to dimer molecules. In this limit, the expression 
obtained for ps reduces to the usual Josephson relation 
for a Bose superfluid in Eq. ([1]). 

As with the derivations for a Bose superfluid in Refs. |6|, 
[rllsllglfiol. the derivation given below for a uniform system 
is based on exact two-fermion propagators, and is not 
approximate. 



II. PRELIMINARIES 

Consider a two-component Fermi gas (e.g., neutral 
Fermi atoms prepared in two different hyperflne states) 
with s-wave interactions between the two components, 
described by the Hamiltonian density (in this paper, fi 
and also the volume V are set to unity) 



V2 

X] I p i^aix) 



2mp 

[/o-0T (a;)'0i (a;)'0T i^) 



(2) 



The two components are denoted by a =t:i and mp is 
the fermion mass. The use of a momentum-independent 
pseudopotential interaction Uq leads to ultraviolet di- 
vergencies that are regularized in the usual way by the 
Lippmann-Schwinger equation. 



1 



mp X ^ mp 



(3) 



Here, is the s-wave scattering length. The entire BCS- 
BEC crossover can be probed by "tuning" the s-wave 
scattering length from small and negative (BCS limit) 
through unitarity (|as| = oo), and finally into the BEC 
limit where is small and positive. 

Although only the case of s-wave interactions between 
fermions is considered in this paper, the analysis given 
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below can be extended to deal with a more general pair- 
ing interaction, as well as Hubbard-type Hamiltonians 
that describe fermions in a lattice. 

The derivation of a Josephson relation for Fermi super- 
fluids given in this paper makes use of the structure of 
the grand canonical thermodynamic potential il.[vs] of a 
current-carrying Fermi superfluid to identify the change 
in the free energy of a superfluid when a velocity is im- 
posed on the condensate order parameter. The superfluid 
density for a superfluid with velocity Vg is then obtained 
from 



(4) 



It can be shown-^ that this is equivalent to the standard 
definition— 



d^F[vs] 



Q2y2 



(5) 



Vs=0 



given in terms of the free energy F — + fin. One can 
also prove (see Appendix A in Ref. [3) that Eq. ([4]) is 
equivalent to = p — p„, where p = mn is the total 
mass density and the normal fluid density pn is given 
in terms of the transverse current correlation functioni^ 
Explicitly, for superfluid flow along the z-axis. 



= p - m{JzJz), 



(6) 



where Jz is the component of the total current operator 
in the z-direction. 

Microscopically, the thermodynamic potential is given 
by the partition function Z, 



-/3^MnZ, 



(7) 



where f3 = (kBT)^^. Functional integral techniques al- 
low us to express the partition function as a functional 
integral over fermionic Grassmann fields and 'ip as^ 



-SliP,i>] 



The imaginary-time action 5 [■(/;, in Eq. 
component Fermi superfluid is given by 



(8) 



for a two- 



(9) 



where the Hamiltonian density H is given by Eq. 12]). 
Here, x = (r, t) is used to denote the spatial coordinate 
r and the imaginary time t = it, and / d'^x = dr J dr. 

A central aspect of the analysis in this paper (and that 
of Josephson) is the existence of a broken-symmetry order 
parameter in the superfluid phase. For Fermi superflu- 
ids, this order parameter Ag is given by the anomalous 
average 



An 



Uo 



E 



(ci,-fcCT,fc)- 



(10) 



Here, c^^k is the Fourier-transform of the Fermi Grass- 
mann field ipcy{x), 



(11) 



where k = (k, Wn) is a 4-vector for the momentum k 
and Fermi Matsubara frequency cij„ = 7r(2n + l)//3, n = 
0, ±1, ±2, and fc • a; = k • r — w^r. 

In order to introduce the bosonic order parameter into 
the partition function Z, the following identity is used: 



I?[A,A*]x 



exp 



{ 



d^x 



Uo 



(A*i/;x7^i + Av;iv;x)]}.( 



(12) 



Substituting this into Eq. ([5]), one obtains the result 
Z = I I?[^,v;p[A,A*] X 

exp { ~ y d^X^'^'tpaix) (^dr ~ ~ 1pa{x) 



Uo 



} 



(13) 



This integral identity (a Hubbard-Stratonovich transfor- 
mation^^) only introduces an auxiliary Bose field A(a;), 
and no approximation has been made. It is straightfor- 
ward to show that the static, uniform component of this 
field gives the order parameter defined in Eq. (|10p . That 
is, using the partition function in Eq. (|13p . one can show 
that 



{A{x)) 



Uo 



^{ci,-kC]^k) = Ao, 



(14) 



where the equilibrium average (A(a;)) is defined by 



(A(x)) 



I?[V',V']2?[A, A*]A(a;)e-'^'=«. (15) 



It is important to emphasize that while the BCS order 
parameter Aq is usually calculated in the mean-field BCS 
approximation, it is not an inherently mean-field quan- 
tity. The "0" subscript on Aq only denotes the fact that 
the order parameter is related to the average occupa- 
tion (macroscopic in the superfluid phase) of a pair state 
with zero total momentum. One obtains the mean-field 
approximation for Aq if the expectation value (• • •) in 
Eq. is evaluated using a mean-field expression for 
the partition function. Here the full partition function is 
used, so Aq is the exact value of the order parameter. 

Having established the relation between the auxiliary 
Bose field A(a;) and the order parameter in Eq. (|14p . 
A(x) can be separated as 



A(x) = Ao + A(a;), 



(16) 
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where A(a;) represents the fluctuations out of the static 
Bose-condensed pair state. 

The partition function given by Eq. (jl3p and the iden- 
tity in Eq. (jl6p will be used below to analyze the super- 
fluid density in a current-carrying Fermi superfluid. First 
we examine the pair fluctuation propagator that describes 
the dynamics of the Bose field A(a:). 



III. THE PAIR FLUCTUATION PROPAGATOR 

The key correlation function of interest in the study 
of the dynamics of Cooper pairs is the 2x2 matrix pair 
fluctuation propagator L(a:, a;') that describes the propa- 
gation of the Bose fleld A{x)^ It is conveniently defined 
in terms of its inverse L~^, with matrix elements given 
by 

(L-i)n(x,x') = -^S{x-x') + Kn{x,x') 

= (L-^h2{x',x) (17) 

and 

{L-%2{x,x') = Ku{x,x'), 
(L-')2i{x,x') = K2i{x,x'). (18) 

Here the matrix two-particle Green's function K is de- 
fined by2i 



K{x,x') 



($(a;)f(x')) ($(x)$(a;')) 
($(a;)$(a;')) (^x)^{x')) 



(19) 



where 



$(a;) = V;t(x)^i(x), ^x) = ij^{x)iP^{x). (20) 



Using Eqs. (flTl) and p8|) . the pair fluctuation propagator 
is given explicitly by 

Uo 



(1 - UoKn){l - U0K22) - UiKi2K2i 

1 - U0K22 U0K12 
U0K21 1 - f/oi^ii 



(21) 



The pair fluctuation propagator can be viewed as the 
propagator for a single composite boson (the Cooper 
pair) and hence, as an analogue of the single-boson 
Green's function D. It should be stressed, however, that 
it is distinct from the two-particle Green's function K, as 
can be seen from Eq. (PT|) . Nevertheless, as discussed in 
Ref. evaluating the pair fluctuation propagator at the 
mean-fleld BCS level (equal to the BCS approximation 
for the many-body T- matrix in Ref. ,25;) , one can show 
that L reduces (up to a constant) to the single-particle 
Bose Green's function [within the Bogoliubov approxi- 
mation; see Eq. (|5D|) ] for a dimer condensate in the BEG 
limit of the BCS-BEC crossover. From this point of view, 
L is seen to be the natural analogue in a Fermi superfluid 



of the Bose Green's function D. Further discussion of the 
relation between the pair fluctuation propagator and the 
two-particle Green's function K can be found in Ref. [23l . 

Motivated by the similarity to the Bose Green's func- 
tion D, we expand the Fourier transform L^^(q, t'm) of 
the inverse pair fluctuation propagator in powers of q as^ 

(L~^)ll(q,Mm) = {L^^)22{<l,-l^m) = Ail^rn- B - 

(L-l)i2(q,z^™) = (L-l)2i(q,i^™) = -i^-Fq2 

+Oiq\>yl), (22) 

where Vm = 27rm//3, m = 0, ±1, ±2, ... denotes the Bose 
Matsubara frequencies. Note that the off-diagonal el- 
ement {Ij^^)i2 has the symmetry (L""'^)i2(q, — i^m) = 
(L~^)i2(q, i/m) and hence the absence of a linear term 
in the Matsubara frequency in its expansion. 

Using the expansion in Eq. ([2^ . the pair fluctuation 
propagator becomes 



L(q,j/„i) = 



1 



+ — + —a^ - 

-D_F 2 ■ 
A a'^ ^'^n 



A 

1 + A + ^q 



where we have deflned the poles of L as 



1 



(23) 
(24) 



This shows that in order for the poles of the pair fluc- 
tuation propagator to be gapless, we must have B = D, 
giving a Bogoliubov- Anderson mode with velocity 



1 



v^j^2BiC-F). 



(25) 



The condition for a gapless mode to exist can also be 
written as 

(L-i)n(0,0) + (L-i)22(0,0)- 

(L-i)i2(0,0)-(L-i)2i(0,0) = 0, (26) 

a result that will be made use of in deriving the Josephson 
relation in Sec. HVl 



IV. THE CURRENT-CARRYING SUPERFLUID 

We now consider the properties of a current-carrying 
superfluid. To introduce a finite superfluid velocity Vg, a 
"phase twist" is applied to the Bose order parameter 
Ao defined in Eq. (fTO)) : 



Ar 



Aoe' 



zm R v = -r 



(27) 



Here, tub = 2mp is the mass of the Gooper pair. Ex- 
panding the Bose pairing field A(a;) about the uniform. 
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static value of the order parameter Aq as in Eq. (|16p and 
applying the phase twist, one finds 



A(a;) Aoe 



vniB'Vs -r 



+A(x) = A(a;) + Ao (e™^"-'- - l) . 

(28) 



Writing the dependence on the superfluid velocity Vs in 
this way emphasizes that the proceeding analysis is based 
on the full Bose pairing field A(x) that includes all fluc- 
tuations about the static order parameter Aq. 

Using Eq. ^E\j in Eq. it is seen that the effect 

of imposing a phase twist on the order parameter is to 
generate a new term in the effective action: 



(29) 



where 



SS 



-Aq I d'*a;[(e™«^=-''- l)$(a;) 

+ (e-™^"=-^- f)$(a;)], (30) 



and $, i" are defined in Eq. (PO)) . Note that the term in 
the effective action Scs in Eq. (fT3|) involving |A(x)p/C/o 
is unchanged by the phase twist to the order parameter 
since the imaginary-time integral over any term linear in 
A{x) vanishes. 

Using Eq. (I29p in the partition function defined in 
Eq. (fT3|). the thermodynamic potential of the current- 
carrying superfluid is found to be 



-S„tt[0]-i5SK 



(31) 



Applying the definition of the superfluid density given by 
Eq. g]) to Eq. ([21]), one finds 



(32) 



Here, (• • •) denotes the equilibrium average in the 
current-free state, given by 



f 



Z[0] 



I I?[V;,^p[A,A*](...)e-^-[ol. (33) 



Note that ({dSS/dvs))\v,=o — 0, by symmetry. 
Using Eq. ^ to evaluate Eq. ^ gives 



Aom 



^ ' d\ (v, •r)2($(x) +$(x)) 



f3 



-Ki2 - K2i]ix,x') 



(34) 



where Vj = Vg/wg and Kij{x,x') denote the elements 
of the matrix two-particle Green's function defined in 
Eq. (fTO|) . Making use of the identity 



($(x)) = ($(x)) 



Ao 



(35) 



Eq. (|34p is naturally given in terms of the inverse pair 
fiuctuation propagator, 



jd'x dV(v,.r)(v,.r')[(L-i)n 



+ (L-i)22 - (L-i)i2 - {\.-%i]{x,x'). (36) 

Fourier transforming this expression, one finds (taking 
Vs = z to lie along the z-axis) 



A 



Ps 



2 q~o dq^ 
-(L-i)2i](q,0), 



lim^ — [(L-^)n + (L-^)22 - (L-^iz 

(37) 



with static matrix elements (L~^)ij (q, — 0). In ar- 
riving at this result, a gapless Bose excitation spectrum 
has been assumed, using the result in Eq. (|26|) . 

Using the expansion in Eq. (I22p to evaluate the second- 
order derivative in Eq. p7p . we see that 



1' 



lim^^[(L-i)n + (L-i)22-(L-i)i2-(L-i)2i](q,0) 



= 4F - 4C. (38) 

This allows us to more compactly write Eq. ([37|) as 

- 2Alm%{C-F). (39) 

Now, from the static (1,1) matrix element Lii(q, 0) in 
Eq. ([23|) . one also finds 



1 [D^-B^ + 2{DF-BC)q^] 

— lim 



lim - , , — xiiii 
q^o q2iii(q,0) q^O 



Bq2 



(40) 



Assuming that L has a gapless excitation spectrum (such 
that B = D), this reduces to 



lim 



1 



2F - 2C. 



q~o q2Lii(q,0) 
Comparing Eqs. ((39)) and (|4T|) . one finally obtains 

Ps 



(41) 



hm 



q^o q2Ln(q,0) 



(42) 



This expression gives the precise analogue for Fermi su- 
perfluids of the Josephson relation for Bose superfluids in 
Eq. We see that the single-particle Green's function 
D for bosons has been replaced by the pair fluctuation 
propagator L, and the square of the BCS order parameter 
Ag plays the role of the square of the order parameter 



<I'2^|(V)P=n, 



(43) 



of a Bose superfluid. 

For a Bose superfluid, Eq. (HS)) gives a simple relation 
between the order parameter $0 and the condensate den- 
sity ric, and these two quantities can be interchanged in 
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the Josephson relation in Eq. ([T]). This is not the case in 
a Fermi superfiuid, however, where the condensate den- 
sity is not a simple function of the order parameter Aq. 
This can be seen from the mean-field expression for the 
condensate density in a BCS superfiuid, given by^^ 

1 

/32 



E 



^mYl Go^2i(k,w„)G'o,i2(k,tj;), (44) 



where Gq is the mean-field 2x2 matrix BCS Green's 
function. Equation ((42|) emphasizes the direct role of the 
order parameter in Josephson's relation, in contrast to 
the indirect role played by the condensate density. 

Equation (I42p gives an exact relation between the su- 
perfiuid density ps and the order parameter Aq in terms 
of the static pair fluctuation propagator L(q, 0). ft can 
immediately be used to study superfluidity in Fermi su- 
perfluids. In Sec. |Vl the Josephson relation is studied 
within the BCS approximation for the pair fluctuation 
propagator. We see how the resulting expression for the 
superfiuid density reduces to the Landau formula for BCS 
quasiparticle excitations. 



V. RELATION TO LANDAU'S FORMULA FOR 
A BCS SUPERFLUID 

An important check of the Josephson relation for Fermi 
superfluids is that it reproduces Landau's well-known for- 
mula for the superfiuid density when the normal fluid is 
comprised of BCS quasiparticle excitations,— 

5/ ,2 



Ps= P 



k 



(45) 



Here, / = [exp(/3i?k) + 1] ^ is the Fermi thermal distri- 
bution for BCS quasiparticles of energy = 
and kz is the z-component of k. The total mass density 
p is given by 



A2 

^0 



k 



l-|(l-2/) 



(46) 



The result given by Eqs. (1^5]) and for the super- 
fluid density is mean-field insofar as it ignores the contri- 
bution to the normal fluid arising from bosonic collective 
modes, as discussed in Refs. [19,25. Consequently, we 
should be able to reproduce Eq. (|45|) by evaluating the 
Josephson relation in Eq. using a mean-field BCS 
approximation. 

Evaluating the pair fluctuation propagator L^^(q, 0) 
within the BCS mean-field approximation amounts to 
evaluating the two-particle Green's function defined in 
Eq. (fTO|) as a loop of two single-particle mean-field BCS 
Green's functions: K = ^ CqGo (schematically). Ex- 
plicitly, Eqs. ini) and (HI]) become 

(L-i)n(q,0) = (L-i)22(q,0) = 

"TT-^ E Go,ii(k,cj„)Go,22(k-q,w„) (47) 
^0 P 



and 

(L-l)i2(q,0) = (L-l)2i(q,0) = 

E G'o,i2(k,tJ„)Goa2(k- q,cj„)- (48) 

k,aJTi 

In this approximation, is equivalent to the nega- 
tive of the inverse pair fluctuation propagator M defined 
in Refs. [l9ll29l . Reference [2^ showed explicitly that the 
poles of L describe the gapless (i.e., B — D) Bogoliubov- 
Anderson spectrum at small q throughout the entire 
BCS-BEC crossover. Furthermore, the combination of 
inverse matrix elements 



(L-i)n(q,0) - (L-i)i2(q,0) = -(C- F)q2 



(49) 



that enters the expression for the superfiuid density in 
Eqs. ([57)1 and (1551) proportional to the static inverse 
propagator for the phase fluctuations of the order pa- 
rameter.^" 

Evaluating the second-order derivative of Eqs. ((Tf|) and 
(|^5)) with respect to g^, after some lengthy but straight- 
forward algebra, one finds 



C-F 



E 



1 

4^ 



1-2/ , df 



2E^ 
^dEl4El 



A. 



dE^ 

2 



3(k) 



where (7(k) is defined by 







2 

) 


r A^l 

1-3^ 








I Eli 



(50) 



(51) 



Note that Eq. ((50)) is the finite-temperature generaliza- 
tion of the Q coefficient defined in Ref. [H- 

Substituting Eq. into Eqs. (gT]) and (021) gives the 
following mean-field expression for the superfiuid density: 



= 2m2 



A2 



1-2/ , df 



2Ek dEj, 



.9(k) 



+2E 



dEl 



El 



(52) 



Using 



df mpEi^ df 



dEk fc^k dk 
and integrating by parts, Eq. ([5^ can be rewritten as 



(53) 



2m 



2 



A 



^ 2E, 



-iV(l-2/).9(k) + 2 



A^ 



df 



V Et dEi, 



2V^ 



1-2/ , df 



2Eu 



A2 
^0 



dE. 



+2™^^E^' 

k 



1-2/) 



3. 



3A| 
El 



,^A^eg df 



k 



Ei 



dEu 



(54) 
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Applying Eq. ((53)) to the last line and integrating by parts 
again, one finds 



A2 



1-2/ , df 
2El 



+ — - 2 



(55) 



The integral in the second line vanishes exactly and our 
expression for the superfluid density reduces to 



1-2/ , df 



2Ei, dEk 



(56) 



Rearranging the mean-field expression for the mass 
density p in Eq. (j46p using integration by parts and 
Eq. ([55)1 (see the related discussion in Ref. jl^, p can 
be written as 



-mp fc; 



(1-2/) 



= 2 



1 - 2/ 



E^dEk 



kl (57) 



Combining this result with Eq. we see that it re- 

duces to Eq. ()45p . Thus, evaluating Josephson's relation 
in Eq. ()42p using the mean-field approximation given by 
Eqs. (|47| and (|48| gives us Landau's formula for the su- 
perfluid density in a BCS superfluid. 

In Sec. IVIl we employ the same mean-field approxi- 
mation used in this section to show that the Josephson 
relation for a Fermi superfluid reduces to the analogous 
expression given by Eq. ^ for a Bose superfluid, in the 
EEC hmit of the BCS-BEC crossover. 



VI. JOSEPHSON'S RELATION IN THE BEC 
LIMIT 

An obvious feature of the Josephson relation for a 
Fermi superfluid is that it must reduce to Eq. ^ in 
the BEC limit of the BCS-BEC crossover, where the 
Cooper pairs are tightly-bound dimer molecules. In this 
limit, where the s-wave scattering length as is small and 
positive, the chemical potential becomes large and neg- 
ative, roughly equal to half the dimer binding energy:'^^ 
p = —l/2mpa1. In this case, |^| ^ Aq, fc^T; / ^ 0, and 
Eq. (|44)l can be solved analytically to give the condensate 
density of dimer molecules, ^^'^^ 



that the inverse pair fluctuation propagator L ^ reduces 
to^^iH 



L ^(q, t^„0 



^)D-i(q,.„), (59) 



where 



D-l(q,M 



^c^mol 



~ ^c^mol 



(60) 



is the inverse single-particle Green's function for the 
Bose-condensed dimer molecules, analogous to the 
Green's function that enters Eq. Here, £q — /2mB 
while C/moi = 4,'!T{2as)/mB is the mean-field interaction 
between dimers, which predicts a dimer scattering length 
of flmoi = 2a^ instead of the exact result a,„oi = 0.6asi^ 
Substituting Eqs. ()55)) and ([55)1 into Eq. (li^ . one imme- 
diately obtains the Josephson relation in Eq. ([T]) for a 
condensate of dimer molecules. 

Of course, one expects Eq. (|^ to reduce to Eq. ^ in 
the BEC limit at any level of approximation, and not just 
at the mean-field level at which Eqs. ([59)1 and ([60)) have 
been derived. The fact that Eq. (|59p is a mean-field result 
only means that the molecular self-energies I]™°'(q, I'm) 
that enter the dimer Green's function in Eq. ([50)1 
are mean-field. Explicitly, writing down the exact single- 
particle Green's function for a dimer molecule, 

D-l(q,z/™) = 



mol 



"il^rn — £q + Mmol — '^11 



"■^12 



^21° ~iVm — £q + A*mol — "^22 



Eq. ([50)1 corresponds to the result 



mol 



and 



mol j 7 

(61) 
(62) 
(63) 

We anticipate that, going past the mean-field approx- 
imation used to obtain the results in Eqs. ([62)1 and ()63p . 
one still arrives at the identity given by Eq. ()59p . except 
that the self-energies will incorporate beyond-mean-field 
contributions. Using Eq. ([CT)) in Eq. ((T)) shows that a 
momentum-independent self-energy will always lead to 
the well-known mean-field result, ps ~ mBnc{T). As first 
shown in Ref. ,25, including contributions from fiuctua- 
tions, the superfluid density in the BEC limit is actually 
given by Landau's formula for a normal fluid comprised 
of gapless Bogoliubov excitations of the BEC of dimer 
molecules. 



k 



A§(r) 



A^(r). 



(58) 



Within the same mean-field approximation [given by 
Eqs. ([37]) and (O], one can show in the BEC limit 



VII. SUMMARY 

Josephson's relation for Bose superfluids gave the first 
explicit identity connecting the two key order parame- 
ters in the theory of superfluids: the broken-symmetry 



7 



order parameter and the superfluid density. It is re- 
markable that such a simple relation exists between two 
such different quantities: the superfluid density that de- 
scribes the response of a system to a transverse current 
probe, as in Eq. ([6]), and the order-parameter in a Bose 
superfluid, associated with the macroscopic occupation 
of a single-particle state. Extending the recent analy- 
sis by Holzmann and Baym^S the analogous identity has 
been derived for a two-component s-wave Fermi super- 
fluid. This gives an exact relation between the superfluid 
density Ps and the BCS order parameter Aq in terms of 
the infrared limit of the static pair fluctuation propagator 
L(q, J^m = 0). 

Using mean-field BCS theory to evaluate the pair fluc- 
tuation propagator, we have seen that the Josephson re- 
lation derived in this paper reduces to the Josephson re- 
lation for a Bose superfluid in the BEC limit and also 
Landau's formula for the superfluid density of a Fermi 
gas with BCS quasiparticle excitations. At first glance, it 
might seem surprising that the Josephson relation-which 
expresses the superfluid density in terms of the propaga- 
tor for collective phase fluctuations-maTia,ges to repro- 
duce this Landau formula for a normal fluid of single- 
particle Fermi BCS excitations. However, this propaga- 



tor is actually a correlation function for the gradient of 
the phase of the order parameter)^ and consequently, 
is directly related to the current correlation functioni^ 
In turn, it is well-known that Landau's formula can be 
obtained by a direct evaluation of the longitudinal and 
transverse components of the current correlation function 
within the BCS approximation.^^ 

The simple structure of the Josephson relation derived 
in this paper should simplify the calculation of the su- 
perfluid density using the standard tools of diagram- 
matic perturbation theory developed for the BCS-BEC 
crossover proble m^^i^^i'^^ to evaluate the pair fluctuation 
propagator L. It also opens the way to giving a rigor- 
ous analysis of the superfluid transition in Fermi systems, 
along the lines of those given for Bose superfluidsfi'^'^ 
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